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A law is designed for simultaneous control of the orientation of an Earth- 
pointing spacecraft, the energy stored by counter-rotating flywheels, and 
the angular momentum of the flywheels and control moment gyroscopes 
used together as an integrated set. of actuators for attitude control. Gen- 
eral, nonlinear equations of motion are presented in vector-dyadic form, 
and used to obtain approximate expressions which are then linearized in 
preparation for design of control laws that include feedback of flywheel 
kinetic energy error as a means of compensating for damping exerted by 
rotor hearings. Two flywheel '‘steering laws 7 ' are developed such that torque 
commanded by an attitude control law is achieved while energy is stored or 
discharged at the required rate. Using the International Space Station as 
an example, numerical simulations are performed to demonstrate control 
about a torque equilibrium attitude, and illustrate the benefits of kinetic 
energy error feedback. 


INTRODUCTION 

Flywheels offer great promise for mincing the mass and extending the lift' of spacecraft: they store 
more energy per unit of mass and last significantly longer than chemical batteries. Moreover, fly- 
wheels can simultaneously store energy and exert torque on a spacecraft, making it possible for one 
system of flywheels to replace two separate systems typically used for energy storage and attitude 
control. When the' mass of tin* two conventional systems is taken into account, the specific energy of 
flywheel systems is expected to he 1 5 to 10 times greater, according to Ref. [1]. lhe attitude control 
system typically represents 11% of the mass of a spacecraft, and batteries make up (> { Z of the* mass, 
therefore a l\\% to 15% reduction in total spacecraft mass is possible by using flywheels. Secondary 
benefits occur as well; since flywheels have higher efficiencies than batteries, a reduction in solar 
array size and mass becomes possible, and reboost propellant can be reduced because the smaller 
arrays product' less drag. Flywheel systems are expected to last 15 years (Ref. [1]) or more whereas 
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typical batteries last only 5 years. The greatest advantage of flywheels over batteries accrues in 
low earth orbit where eclipse happens more frequently and for longer periods than in higher orbits: 
repeated charge and discharge cycles, and high depth of discharge significantly degrade batteries 
over time. 

Consideration was recently given (Ref. [2]) to replacing the International Space Station (ISS) 
batteries with a Flywheel Energy Storage System (FESS): hence, there naturally arose the thought 
of using the flywheels to assist the ISS Control Moment Gyroscopes (CMGs) in controlling attitude. 
A numerical investigation of the merits of this idea requires a feedback control law designed for 
CMGs and flywheels used together as an integrated set of effectors. The current CMG control law 
described in Refs. [3] and [4] minimizes a cost function involving spacecraft attitude and angular 
speed, and CMG angular momentum. A CMG steering law such as the one developed by Hans 
Kennel determines the speeds of the two gimbals (in which each constant speed CMG rotor is 
mounted) needed to produce the torque requested by the control law. As a natural extension of 
the present approach, we seek a new control law derived from a cost function that includes flywheel 
angular momentum in addition to the aforementioned quantities. Also needed is a 4k fiywheel steering 
law," a counterpart to the CMG steering law that will determine the motor/geuerator torque to bo 
applied to eac h member of a counter-rotating flywheel pair such that rotational kinetic energy is 
stored or discharged in the required manner, and the net torque 1 requested by the control law is 
produced simultaneously. 

A review of the literature does not reveal any suitable existing control laws. In Ref. [h], Notti. 
Cormack. and Klein give a sketch of a control law and an energy distribution law: however, this work 
is not applicable primarily because each flywheel rotor is assumed to be supported by two gimbals 
and the ISS FESS did not contain any gimbals. In addition, the control law lac ks flywheel angular 
momentum as a feedback parameter. Recent work on control laws for integrated power and attitude 1 
control systems deals either with sets of four or more 1 flywheels whose 1 spin axes arc 1 uon-collinear 
(Refs, [(if [8]). or with sets of variable-speed, single-gimbal control moment gyroscopes (Refs. [9] 
[11]): neither of these type's of configurations are directly applicable to the 1 counter-rotating flywheel 
pair arrangement of the FESS. Varatharajoo and Fasouias [12] develop control laws for a spacecraft 
using a counter-rotating flywheel configuration: however, they only consider a pitch-axis controller. 

Hall’s control law. proposed in Ref. [(]]. is an open-loop scheme (rather than a feedback law) for 
performing large-angle attitude maneuvers. It does not account for gravitational and aerodynamic 
torques which have a significant effect on the motion of ISS, and therefore can not be used for 
maintaining torque equilibrium attitude, the primary job of the CMGs. Hall introduces a flywheel 
steering law based on a matrix pseudo-inverse 1 ; it is applied in each of Refs. [7] [11]. Tsiotras. 
Shell, and Hall employ Lyapunov stability theory in Ref. [7] to develop a ft**! back control law which 
performs well in simulations involving disturbance 1 torques: however, flywheel momentum is managed 
by expenditure of propellant. Costic et al. develop in Ref. [8] a nonlinear controller which includes 
an adaptive scheme for estimating the 1 mass distribution of the spacecraft, but they do not address 
momentum management. 

In Ref. [9]. Fausz and Richie extend the work of Hall to a nonlinear feedback controller applicable 
to a set of variable-speed, siiigle-giinbal control moment gyroscope's. Together with Tsiotras. they 
continue their discussion in Ref. [10] and present simulation results, but momentum management is 
not addressed in either of the two papers. Yoon and Tsiotras develop an adaptive nonlinear c ontrol 
law in Ref. [11]. and incorporate wheel-speed equalization to reduce the possibility of singularities 
and keep the wheel speeds within acceptable limits. Numerical simulation results show that attitude 1 
and power profiles can be tracked even when the spacecraft, inertia properties are unknown. Of all 
tlu 1 papers mentioned heretofore. Ref. [11] is the only one in which attitude control, momentum 
managememt. and power management for flywheels is considered in an integrated fashion. 

It is important to note 1 that Re'fs. [Ci] - [12] fail to take 1 into account damping torque exerteM by 
the 1 spacecraft and a ftywhcvl rotor on (‘ach other: in practice, this will cause 1 the 1 actual rotational 
kinetic energy possessed by the 1 flywheels to differ from the 1 requires! amount. None 1 of those works 
include 1 any feedback of errors in flywheel power or kinetic* energy, somi'thing which must be 1 done 1 
under realistic conditions. 
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Figure 1: Spacecraft with Flywheel and CMC 


In what follows we present equations of motion to be used in numerical simulations and in control 
law design, develop two flywheel steering laws, design an algorithm for managing momentum and 
maintaining torque equilibrium attitude, and present, simulation results showing the performance of 
the control law as well as the benefits of flywheel kinetic energy error feedbac k. 

EQUATIONS OF MOTION FOR SPACECRAFT WITH FLYWHEELS 
AND CMGS 

Dynamical Equations 

The system of interest, 5, is composed of a rigid body B moving in an inertial or Newtonian reference 
frame' N, and several rigid axisymmetric rotors /?].... <R P whose mass centers are fixed in B. A 
subset of the rotors , Rjr have spin axes fixed in B so that these rotors represent non-gimballed 

flywheels or reac tion wheels. Each of the remaining rotors R f) are attached to B with one 

or more massless gimbals which permit the direction of the spin axis to change relative* to B: these 
rotors thus represent a number of CMGs. C = p - T. (The latter subset could contain gimballed 
flywheel rotors as well as CMGs, but we concern ourselves in this work only with non-gimballed 
flywheels.) This system is illustrated in Fig. 1. with rotors /? L >. ■ • • omitted for the sake* of 

clarity. 

The equations of motion are derived using Kane's equations (Eqs. (0.1.2), Red. [13] ) 

F r + F* =0 (r = L... ,«) (1) 

where 1 F, are generalized active force's for S in A r . F* are* generalized inertia forces for S in A , and 
n is the number of degrees of freedom of S in N. 

The system S is holonomic and therefore' a complete description of the* motion of -S’ in N requires 
ii generalized speeds mj u n , conveniently chosen as follows. 

* uj n — u i b j + U )h ) + (2) 

where A uj H is the angular velocity of I? in A\ and bj, bj. and b* an* a set of mutually orthogonal, 
right-handed unit vectors fixed in B. 



Wo introduce unit vectors (3 t fixed in B such that they are each parallel to the spin axis of 
a flywheel rotor 77, , and therefore to the angular velocity Rl of R } in B. Generalized speeds 
u.i. ■ • • • ujF+a associated with the flywheels can then be used to write the angular velocities as 

V' = «, +;i / 3, (i-- i T) (:i) 

The inner giinbal of each CMG is fastened to B with a revolute joint in a single-gimbal config- 
uration. whereas a double-gimbal CMG has the inner giinbal attached with a revolute joint to an 
outer giinbal. which is in turn mounted in B with a second revolute joint. The axis of each revolute 
joint is assumed to pass through the mass center of the rotor, which is thus fixed in B . In order to 
maintain generality with regard to the number of CMG rotors and gimbals, the angular velocities 

Rf * 1 B (j 0 H( ‘ of the rotors relative to B are not written explicitly, but they must be functions 

of the generalized speeds t/^+ 4 , . . . , 3 , where one generalized speed is required for every gimbal. 

The final three generalized speeds are associated with the velocitv of the mass center S* of S in 

A\ 


A V s = ■u„_2 h] 4- U fi - in 2 4- u n (4) 

when* no, and h;j are a set of mutually orthogonal, right-handed unit vectors fixed in N. 

Let a be the set of forces exerted on S except those* exerted by B and R } on each oth(T (/ = 

1 p). The forces in a acting on B,R 1 R p are equivalent to single forces F/^.Fi F /f 

applied at the mass centers B\R\.... .77* of bodies B,R j. .. .R f) . respectively, together with 
couples whose torques are M/*, M| M ; ,. 

To account for the forces exerted by B on R { . we regard them as equivalent to a single force* 
F R ' r ' applied at R * , together with a couple whose torque is WL B/ R ‘ . Sinc e 77* is fixed in B. 
contributes nothing to F, (i — 1 p:r = 1 , . . . , n). 

The* generalized active forces for S in N are obtained by application of Eq. (4.(j.l) of Ref. [13], 
and are given by 


F r = N v?' • F + -M + J2 • (m, + M ,,/w ' ) (r = 1 n) (5) 

?= 1 

where A v* is known as the rth partial velocity of S * in TV, is the rth partial angular 

velocity of B in N. and so forth. The vector F is the resultant of the forces in a acting on S. 
F — F/? 4~ Yli=\ F/- and M is the moment of a about $* . given by M = M# 4- r lS * D * x F/j 4 - 
! (M, 4- r s R ‘ x F,). where* r lS B is the position vector from S * to B* . etc. 

The generalized inertia forces for S in N are formed according to Eqs. (4.11.5) (4.11.7) of 
Ref. [14], and are given by 
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where ni s is the mass of system S. A a tS is the* acceleration of S * in A\ I lS,/ ‘ s is the inertia dyadic 
of S for S\ A 'a /} is the angular acceleration of B in A r , \ R,/R > is the inertia dyadic of R, for 77*. 
and the* central angular momenta of 77, in B . and in A r , are denoted respeetivelv bv R H R< ? R > . and 
A H R /R * . Difl’c TcMitiation with respect to time in B. and in N. arc* indicated respectively by B d/dt 
and A d/dt . 
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Ac cording to Eqs. (1) the generalized inertia forces from Eqs. (G) are added to the generalized 
active forces from Eqs. (5) to yield veer tor-dyadic equations of motion for a spacecraft containing 
flywheels and CMGs 


A 'v; s • (f - ms A a ,s ) 


+ { M - 


r sy.s- . N a H + N u B xI s/s- . JV W n 


\ L 

+ Z f X b h w ‘cA + ST (V x !«•/«: X «u» w ') 

,=i V / »=^+i v 


v u> /J 




+ ]T "w/' 1 • | M> + M H/R ‘ - — ) - 0 (r = 1, . . . ,«) 


/= i 


Eqs. (7) are completely general with regard to the number and orientation of flywheel rotors, and 
the number of CMC rotors and gimbals. In this form, they are applicable to variable-speed CMGs. 
These equations of motion, and the expression for generalized inertia forces, can be compared briefly 
to previous work. 

Reference [14] is concerned with gyrostats and relevant equations that can be 1 dealt with easily 
by an analyst, and quickly by a computer. Expressions for generalized inertia forces are presented 
separately for a gyrostat containing a single cylindrical rotor, and for one containing a single' spherical 
rotor; an underlying general relationship (CG5) developed in Appendix C can lx* shown to give rise 
to Eqs. (G) presented here, when no CMGs are present (C = 0. thus p — JT). The term in the 
second line* of Eqs. (COG) accounts for but a single rotor, although additional rotors can be handled 
straightforwardly by adding a sum of similar terms. The correspondence between Eqs. (CG5) and 
our (G) is shown by appealing to Eqs. (24), (CGI), (C35), (C24), (4), and (C34) of Ref. [14], and 
replacing their labels G. D. and A for the gyrostat, rotor, and carrier respectively with our S, /?,. 
and D. The use of the system mass and inertia scalars, together with the moment of inertia for the 
axis of symmetry of each rotor, is shown in Ref. [14] to lead to greater efficiency than use 1 of mass 
properties of individual bodies in a gyrostat; this advantage happens to accrue to Eqs. (G) and (7) 
developed hen*. 

Rheinfurth and Carroll present in Ref. [15] vector-dyadic equations of motion (27) for a spacecraft- 
composed of a rigid carrier and a rigid appendage whose mass center is fixed in the carrier. Additional 
appendages are accounted for easily by forming a sum, as they do in Ecj. (l(i), but the resulting vector- 
dyadic expression will give rise to only three scalar equations. It is pointed out near the bottom of 
p. G that a CMG can lx* regarded as an appendage; however, the motion of every appendage relative 
to the carrier must lx* prescribed if the 4 equations are to serve as dynamical equations governing 
the motion of the carrier. Reference [15] does not contain counterparts to the n — G of our Eqs. (7) 
that govern the motion of the rotors, or to the 3 equations governing the translational motion of 
the system. It can be shown rather easily that, Rheinfurth and Carroll's Eqs. (27) give* way to tlx* 
first three* of Eqs. (7) here when all rotors are permitted to be* CMGs (JF = 0. thus p = C). After 
forming the required sum, and replacing them* symbols L with M, \ with I s /S . I p with I /v R ‘ . 
U with N u ! *' u; p with b uj r ' , (U) v with N a B . and (i p ) p with n a R '. it becomes evident that 
forming dot products with the* resulting expression and three vectors A u> R produce's the* first thre*e 
scalar relationships given by Eqs. (7). 


Dynamical Equations for A Complex Gyrostat 

A space'craft known as a simple gyrostat is described in Se*c. 3.G of Ref. [Hi]; the* system S in the* 
preceding discussion becomes a simple gyrostat when the number of fly w heeds T is equal to 1, and 
when no CMGs are* present (C — 0, p — T). A spacecraft with more than one* fiywhe*e*l. such as the* 
one* shown in Fig. 2. will be* referred te) as a ce>mple*x gyrostat; expiations of motion with T — G are 
give’ll in the* following material. 
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Figure 2: Spacecraft with Flywheels 


Without a great loss of generality one can at this point work with six flywheel rotors i?j . . . . . 

[T — (>) arranged in three counter-rotating pairs as shown in Fig. 2. with the spin axes of R\ and 
R.\ parallel to bi, R > and R-, parallel to b 2, and R-j and R parallel to b 4. The generalized speeds 
i/4, . . . , w.y associated with the flywheels are then used to form the angular velocities of R; in 

R (i = 1 (>) 

n u> hl = a.jbi. h uj Rl — Ui'yh), 

h uj Hi = Mr>bi, B u) Hr = i/yb-j. b l o Hu — (8) 

Up to this point the moment about R* of forces exerted by B on R, has been represented by 
M B/ r ; henceforth we regard the dot product M B/ ^' • (3 t (i — 1, . . . , JF) as the sum of two contribu- 
tions. The first is from a motor /generator, and will be denoted by M W//f ' * for convenience. The 
second is due to damping, related by a constant of proportionality C\i to the angular speed of R } rel- 
ative to B. For example, in connection with rotor R \ , M . h/ • b } is replaced by WL B/Ba • b] — (7,/Wo. 
in Ecjs. (7). 

Unit vectors bi, b_>, and b 2 are taken to be parallel to central principal axes of inertia of S , so 
that 

I ^ — /jb]b] + / 2 b 2 b 2 4- / 2 b 2 b 2 (9) 

wliere /[, / 2 , and 1$ are central principal moments of inertia of S. Equations (7) then yield twelve 1 
scalar relationships; upon decoupling the first three 1 from the 4th through 9th. they can be written 
as 

(7i - 2,7)// 1 = (1> - 7:j )"■->".! - .7 [m 2 (mk + w-<i) - </.-«(«« + //?)) + M\ 

- (M, + + M, + • b, + C/(«i + ur , ) (10) 

(I‘2 - 2J)ii'2 = (hi - + «;,) - in(u H + My)] + 

- (m_, + M u/W -' + Mr, -f M /,/Wr -) • b_> + Crf(w (i + u 7 ) (11) 

(h ~ 2 ./)(/,( - (/ 1 - I>)ii\u> - J[(/i («g + (1~) ~ u>(u.i + Mr,)] + M a 

- (m : , + + M,i + M" /w ") • b :i + C’ ( ,(//* + «.,) (12) 
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j(u \ + u.\ ) = 

(m, +M b/K| ) 

• bi - C rf «4 

(13) 

./( // ] -f ilr t ) — 

(m.i + 

• bi - C,i Ur, 

(14) 

J ( h ) T in i) — 

(M , + 

• bo - C,iu r> 

(15) 

.7(7/2 + th) — 

(m 5 + 

• bo - C,IU T 

(16) 

.7 (ii 3 4- //*) = 

(m.j + M /,//r ‘) 

• b;j - C'.ill K 

(17) 

J{*h + '«!)) = 

(m« + M B/W,i ) 

• b;( — 

(18) 

ni$u r = 

F r -n ( r = 10. 

11. 12) 

(10) 


when 1 J is the central principal moment of inertia of a flywheel rotor for its axis of symmetry, the 
scalars M } in Eqs. ( 10)— (12) are defined bv the relationships Mj = M • b j. and Fj in Eqs. (19) are 
defined as F, = F • n, (j = 1.2,3). 

Eqs. (10) (If)) are thus a complete set of nonlinear dynamical equations of motion for a complex 
gyrostat composed of a body D and six axisymmetrie rotors R j,... , /?(> whose mass centers and 
spin axes are fixed in D\ the rotors are arranged in pairs, with the spin axes in each pair parallel to 
each other and to a central principal axis of inertia of the gyrostat. 

If a rotor R, does not possess a magnetic dipole moment, and is housed inside D where it is 
protected from tlu v action of aerodynamic forces, then the principal contribution to the 4 external 
torque M, is gravitational moment exerted by the celestial body about which the gyrostat orbits. 
An expression for an often used approximation of gravitational moment is given in Eq. (2.9.8) of 
Ref. [19]. where it becomes clear that the dot products M> • b, and * b, (z = L2.il) all vanish 

because* in each case the* unit vector b, is parallel to an axis of symmetry of the rotor. \\ ith C\i = 0. 
Eqs. (10) (12) and (17) arc 1 easily shown to reduce to Eqs. (3.7.28) (3.7.31) in Ref. [19] for a simple 
gyrostat in which the spin axis of the single rotor is parallel to b 

Approximate, Linear Equations for a Spacecraft with Flywheels and CMGs 

In Ref. [3], Wie et al. develop a scheme for controlling a spacecraft’s attitude 1 and managing the* 
angular momentum of a collection of CMGs, and the results arc 1 applied to a space station. In the 1 
case* of a spacecraft, carrying CMGs and no flywheels (J- — 0, p = C), the first three of Eqs. (7) can 
be shown to give 1 rise 1 to the six relationships employed as a basis for control law design in Red. [3], 
namely Eqs. (1) and (3) therein. An essential step in the demonstration consists of neglecting 
the 1 second sum (which receives contributions only from CMGs) in comparison to the first sum 
appearing in Eqs. (7), bast'd on the assumption that the 1 CMG gimbal speeds are much less than 
the 1 rotor spin speed. In addition, central moments and products of inertia of S are regarded as 
constant in Red’. [3], based on the assumption that reorientation of CMG rotors (and gimbals) does 
not significantly redistribute system mass. Both assumptions arc 1 quite 1 reasonable in the* case 1 of the 
International Space Station and its four CMGs with constant rotor speeds of 9.900 rpm. 

There are n dynamical equations of motion (7). ti —3 of which govern rotational motions of D and 
the CMG rotors. The first three of these are approximated by the three 1 Eqs. (1) of Red. |3]. together 
with the three additional relationships introduced in Eqs. (3) of Red. [3 j to represent the 1 effects of 
CMGs. This approach doe's not account entirely for all n - 9 of Eqs. (7) that must follow from 
the third sum therein, where v — 9 is equal to the number of gimbals (involute joints) supporting 
the 1 CMG rotors in D ; however, for the purposes of designing a control law the approach c an be 1 
extended to represent the* effects of flywheels in a similar manner. There 1 results nine 1 approximate 1 
dynamical equations for a spacecraft, with fiywheeds and CMGs. To describe the 4 orientation of B 
in a local vertical, local horizontal reference frame* I. one 4 may choose 4 a body-three. 2-3-1 rotation 
sequence as set forth on p. 423 e>f Red. [19]. The 4 sequence is also known as pitch-yaw-roll, with the 1 
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angitis demoted by 0j. Oo. and 0 :i respectively, and there are three associated kineniaticai equations 
of motion. The set of twelve equations can be linearized about an Earth-pointing motion associated 
with a circular orbit of rate n, and written as 


hui — (/a - l 2 )( 3 n 2 0 :i + mi 3) — T\ - T\ + w\ 

(20) 

I 2 Ft 2 — {h ~ /i)( 3 ?r^i) - To - fo 4- wo 

(21) 

= ( h ~ h) 

ml j - 7-3 - 73 -f 

(22) 

hi = nh-.i + t 1 . 

ho = To, I13 — —nh 1 -f Ti 

(23) 

Hy = nH;i 4 - ?\. 

Ho — fo. H;\ = -nH\ + fa 

(24) 

0 , = u 2 , 62 = 

M 

+ 

IIS 

II 

1 

II 

(25) 


where all quantities with a tilde over them are referred to as perturbations, and considered "small.' ” 
The moment M is regarded as the sum of two terms, the gravitational moment exerted on S and 

reflected in the terms with a coefficient of 3n 2 . as well as all other contributions w. with tv, == w * b T 
(r = 1.2,11). The scalars r r are defined in terms of the contributions from CMGs to the first sum in 
Eqs. (7), 


r, £ I J + ,v w B x b h «./«; j . b 

r /?,//> \ / p 

]T b H h,/r * j + B X I ]T /j H R,; i{ 

) \i=;F+ 1 


d 

dt 


•b, (7- = 1.2,3) (2(i) 


The scalars h, and h r denote' the perturbations of the quantities 


h r = h • b, = H ‘ 

\/' = .F+ 1 


fir 


•b r . h r = T h r (r — 1, 2, 3) 
dt 


27) 


The scalars r v . H r , and H } . and the vectors r and H. associated with flywheels, arc involved in 
similar relationships with the limits of the sums going from 1 to T , rather than from T + 1 to p. 

Equations (20) -(25) are supplemented with differential equations governing ( h r dt, f H, dt. and 
/ dj dt (r = 1.2.3) such as those introduced in Ref. [3] to eliminate biases in the time histories of 
the integrands. These relationships become better suited for control law design after they have been 
made nondimensional by use of the following definitions. 





r * A Hi_ 

' ~ h n 


(/’ - 1,2.3) 


(28) 


hn 1 ' 


r, = 


/, »- ' 


h »- 


(i = 1.2.3) 


(29) 


t* = nt (30) 

Different subsets of the 21 uoitdimeusional, linear differential equations can lie written in t lie form 
required for application of the Linear Quadratic Regulator technique. 

{x} = [.4]{.r} + [5]{r} + {ir} (31) 

to obtain various linear controller designs with the* goal of determining the values of r and r that 
best control the orientation of B in L . and minimize the 4 magnitudes of h and H. The column matrix 
{ 7~ } is dimensioned (ix 1 with the* elements 

W = [ T\ T-2 T;{ f,* n t* ]' (32) 

where the* superscript T indicates the transpose of a matrix. 


8 



FLYWHEEL STEERING LAWS 


The great benefit of utilizing flywheels is that they can serve simultaneously as attitude control 
actuators and as energy storage devices; this dual role requires that r. obtained on the basis of 
attitude control considerations, bo applied in a way that allows energy to be stored or discharged 
as needed. A flywheel rotor R, is suspended in a vacuum housing in D with magnetic hearings, and 
relative motion between R, and B is brought about by a motor-generator that enables B to exert 
on R, a torque with magnitude! M h /H) * (3 h the purpose of which is to produce and change rotor 
momentum in order to furnish attitude 1 control, and to alter the rotor s rotational kinetic energy. 
An important measure of energy storage is power, or the rate at which rotational kinetic energy 
is changed. In this section relationships for • (3 t as functions of power, f, and rotor speeds 

are developed: they are referred to collectively as a flywheel steering law because they are similar in 
nature to a CMC steering law that determines gimbal speeds (and thus, indirectly, gimbal motor 
torques) needed to produc e r as requested by a control law. Two such steering laws are presented: the 
first is the* result of simply prescribing the total power of the flywheel system, whereas specification 
of the power required for each of three flywheel pairs gives rise 1 to the second law. Bearing friction 
and damping an 4 neglected in the design of the steering laws. 

The flywheel rotors arc* arranged in counter-rotating pairs; each pair, denoted by F). consists of 
rotors Rj and (/ — 1,2,3). Referring to Eq. (5) of Ref. [17]. the power H P h * of F { in B can 

be* expressed as 

11 P = (l«- *' + (l Wi J • 'V* «>' ■ 1 ) • H u> ' (i = 1. 2. 3) (33) 

or. in view of Ec}s. (S), 

B P^ — J ( U >1+2 u 27 + 2 + h 2 ? + 3 ( 1 2 ? + 3 ) ( ? ~ 1 , 2 . 3 ) ( 34 ) 

and therefore the* total power may be expressed as 

3 3 

Hph ^ ^ Bpl , _ jy (u, 2 ; +2 U2 } + 2 + ^2/ +3 ^2/ +3 ) (35) 

?:= l 7=1 

Now, if ii ) , U‘ 2 , and 1/3 are assumed to be small in comparison to c/.j //<)• Crf is neglected, and 

M, • b, and • b, (i = 1,2,3) vanish for the reasons put forth earlier, then Eqs. (13) (18) may 

be* approximated as 

Jh 2irl ft M H/Hi • b„ * b, (i = 1.2.3) (30) 

which, together with the* counterparts to Eqs. (20) and (27). load to 

(r A u; H x H) • b, - H, — J(ii 2 i +2 + thi+s) 

ss (m R/H ' + M" /w ' ,; ') • b, (/• = 1.2.3) (37) 

or 

M fi/R ‘-' • b, ss (f - N u> ri x H - M w/ "') • b, (i = 1.2.3) (38) 

Substitution from Eqs. (30) into Eq. (35) produces 

11 P ,r « Y. («*+ + U2,+:M H/Hi * ’) • b, (3!>) 

1= J 

Eqs. (38) and (3D) const itiltc a system of four equations. linear in the six unknowns M /? h • b, 
and • b,. (•/ = 1.2.3). This uudordetomiincd system can be written in matrix form as 

[-4,! {«} = {;} (40) 
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where* 



1 

1 0 0 

0 0 




■ . 1 A 

0 

0 1 1 

0 0 




A«\ = 

0 

0 0 0 

1 I 



(41) 


_ U. ! i 

[ir, no iij 

U$ My 




M = [ 

M b/I{ ' • b, M l,/Wj • b| • b > 

M B/Rr ' • b_, ■ b, M B/H " • b. 

] T 







(42) 

{-'} = [ 

(r- ‘ 

v u > B x H) 

• b, (f- 

N u> B x H) 

• b_, (r - A 'w B x H) • b :J B P F ] T 

(43) 


One may solve Eqs. (40) by forming a matrix pseudo- inverse such as the one presented in Ref. [18] 
and developed by Moore and Penrose for under determined systems. 

[A,] + = [As] T ([-4 s ][A,] T ) _1 (44) 

which yields the solution 

{y} = [A H ] + {z} (45) 

that minimizes the sum of the squares of the unknowns. {//} 1 {?/}. (There exist several other 
performance measures that could be considered in solving an underdetermined system of equations.) 
This is in essence the suggestion made by Hall in Sec. 4 of Ref. [0], yielding a steering law wherein the 
power and attitude control requirements are met simultaneously, and a function of the instantaneous 
motor torques, the sum [(M ^ y H ' • b*) 2 + (M^ ' • b,) 2 ]. is minimized. Upon making tin* 

definitions 

(l[ = f /5 — f/ 4 , do = Uj — U(i . d;i = Hi) - (40) 

A A A 

.s i — u-r, 4- «i. s> = a 7 -f Uo, > 4 - Us ( 4 / ) 

the pseudo- inverse can be written explicitly as 

[A,] + = (4«) 



2ur } d[ 

+ dr + d , 3 - 

sod 1 

s^d) 

-2d, 


— 2u±d\ 

+ dr + d,, 2 

—sod\ 

-s.id 1 

2d, 

1 


•S' 1^2 

2ujd> + d\ 2 4- d :i 2 

s :i d) 

- 2d 2 

2 (d\ 2 4 - d> 2 4- d A 2 ) 


— .s [do 

~2uc } d'2 4- d\ 2 4- <7;j 2 

— S:\d> 

2d> 



»i<k 

Sod.i 

2th)d A + d\“ T df 

-2d,, 



s \di 

— S‘)d‘A 

-2 Uxdn 4- c/i “ + do 2 

2d,, 


It is worth noting that if the rotor speed differences c/ 2 , and r/ : * vanish for all three flywheel 
pairs, the pseudo-inverse becomes infinite and the steering law does not furnish a result. Since the 
rotor speeds of the flywheels in each pair will normally have opposite signs, this condition should be 
unlikely. In addition, examination of Eqs. (42), (43), (45). and (48) reveals that when t I k* flywheels 
arc* not required to provide attitude control (r = 0), and the condition of counter-rotation is present 
thus H = 0). this steering law dictate's that M /j + ’ * b, = * b,. and 

counter-rotation is preserved. 

As an alternative to dealing with the four Eqs. (38) and (30). one could replace the single Eq. (:«») 
with three 1 others by choosing to divide the power requirement evenly among the three flywheel 
pairs, B P Fl -- M pi — ± B p r \ yielding six equations in six unknowns. Substitution 

from Eqs. (30) into Eqs. (34) gives the three new equations 

n P F ' « («2i+2M /,/Wi + « 2H; ,M w/w '- : ') • b, (/ = 1.2. .4) (40) 
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Now, substitution from Eqs. (38) gives 

B P K « m 2 , +2 (m b/i< ' ■ b,-) + « 2<+3 (f - N u B x H M b/w ') • b, 

= (« 2,+2 - «»+*) (m h/r ‘ • b,) + m 2 ,-+3 (r - A "u> w X H) • b, (/ = 1 . 2. 3) 
which can Ik* rearranged to yield 


M /i/W| • b, 


H P F ‘ - «2,+.( (T - A ' 


“ A r . . Zi 


up x 


H) • b> 


('' = 1.2.3) 


U2>+2 — ' U 2i+:J 

and. when oik* substitutes from this expression into Eq. (38), tin* result is 

"■>,+■> (f - N u B x H) • b, - 11 P K 
M li/ "• ■ • b, = -^±zl - (» = 1. 2. 3) 

«2?+2 - U,2i+;* 


( 50 ) 


(51) 


(52) 


Ecjs. (51) and (52) constitute another flywheel steering law, indicating the moment that must Ik* 
applied by a motor-generator to each of two rotors belonging to a counter- rotating pair in order to 
apply r as called for by an attitude control law. and at the same time satisfy power requirements 
specified by B P h '. Each flywheel pair is expected to operate with the sign of the rotor speed 
opposite the sign of u?i+ 2 - hence, the denominators “ il 2 i+A should remain well away from 

zero. In the event that the flywheels are not required to participate in attitude control, and the 
condition of counter-rotation is present, the steering law yields M 7>> ' ' * b, = — M /:? !i> * b, and 

thus maintains the* condition of counter-rotation. 


LINEAR CONTROLLERS 

Laws for controlling the orientation of B in L, and the momentum of flywheels and CMGs. have 
been designed using the infinite-horizon Linear Quadratic Regulator technique, in which a scalar 
quadratic performance index given by 

V = P ({./’} ' [Q] {./■} + {r} T [/?] {r}) dt (53) 

is minimiz(*d subject to the* linear Eqs. (31) with {IV } = {()}. The technique yields a state feedback 
gain matrix [K] which in turn is used to obtain {r}. 

{r} = - [A'] {*} (54) 

As sugg(*st(*d by Bryson and Ho in Ref. [If)], the weighting matrix [Q] can be chosen as diagonal, 
and unity should he approximately equal to the product of Qjj and the square of the maximum 
acceptable value* of the associated element .r y of {.r}. Maximum acceptable values of several param- 
eters used in constructing [Q\ an* listed in Table 1. Before constructing \Q}. the values in Table 1 
must be made* nondiniensioiial to be in correspondence with tlu* elements of {./’}■ Likewise, a diago- 
nal form of [7?] is convenient, with unity approximately equal to the product of /?u and the square 
of the maximum acceptable value of the associated element of {r} (k = 1,... .6). The maximum 
expected value of each element of {r} is taken to he 2.7 N-m, which must be nondimensionalized 
before constructing [/?]. 


Simulation Parameters 

Other values required for a numerical simulation are as follows. 

Moments and products of inertia of S with respect to 5* an* taken to be 



50.28 

-0.3!) 

-0.24 “ 

I s/s ' } = 1.3(> x 

—0.30 

10.80 

0.10 

J 

-0.24 

0.10 

58.57 
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Table 1: Maximum Acceptable* Values 


Parameter 
(i = 1,2,3) 

Max Value 

e, 

1 deg 

U t 

0.2 deg/s 

hi 

0.779 N-m-s 

Hi 

0.779 N-m-s 

f h > dt 

2.7 x 10° N-m-s- 

J Hjdt 

2.7 x 10° N-m-s* 

J0,dt 

1000 deg-s 


which are the same values (in metric units) as those associated with Phase 1 in Table 1 of Ref. [3j. 
with the exception of I\j which is suspected to be a typographical error, the correct value being 
-0.24 x lif slug-ft 2 . 

The value of w used in Ref. [3] represents the moment about S * of aerodynamic force's exerted 
on the Phase 1 configuration of S when the attitude is near TEA; it is given by 

w = 1.30[(1 -4 sin nt + - sin2//f)b| t(4 + 2 sin nt + i sin2n/)bj + (1 4- sin nt + ^ sin 2//^ ) b ^ N-ni 

(50) 


where //. the magnitude’ of ;N u > lj . is taken to be 0.001131 rad/s. 

Each Awheel pair in the FESS is required to discharge 4.400 W of power during the portion of 
the orbit that lies within the Earth's shadow, known as the period of eclipse. The remaining portion 
of the orbit, during which sunlight readies the spacecraft, is taken to be twice as long as the eclipse. 
Therefore, for each orbit, the total power that must be supplied by the 48 pairs of flywheels in the 
physical system is given by 

BpF _ f 48 x 2, 200W = 105.0 kW 0 < t < (charge) 

1 ~ \ 48 x -4,400 W = -211.2 kW < t < Izl (discharge) 

where the bar over P indicates a known function of t to be used in connection with the pseudo- 
inverse steering law developed previously. The alternative steering law is referred to as a divided 
I jower law because the power requirement is divided into three equal parts, therefore B P h ' — 
H P F * = B P F " = y j P F .or 


( 35.2 kW ()<t <ri2L (charge) 

\ -70.4 kW < t < "tt* (discharge) 


(' - 1.2.3) 


(58) 


The FESS would have been made up of 90 flywheels; since the present model involves only 0 
rotors, a scaling factor of 90/0 = 10 is used: thus J — 10 x 0.3010 kg-nr — 4.82 kg-nr. 


Energy Feedback 

The power possessed by the flywheels, given in Eq. (35), can differ from the value required by 
Eq. (57) when C\i / o. If left unchecked, such unwanted resistance or damping will lead to a 
difference between B P h and B P h that increases with time, to rotor speeds that exceed their 
maximum and minimum limits, and to singularities in the steering laws. These deleterious effects 
of damping can be eliminated by the* control system through feedbac k of rotational kinetic* energy 
error. 
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The total rotational kinetic energy B K r of the flywheel rotors relative' to D can he expresses 


as 


■A- 4 XX' 


m 


i=\ 


and the power of F in D is given by the derivative of 15 K *' with respect to L 


BpF _ d_B K /■ 


(GO) 


One can regard the 1 required power as the time derivative of a required kinetic energy of F in 

D 


BpF _ Bj-'F 

dt 


( 01 ) 


and define* a kinetic energy error c/. as the quantity 


a- = B K F - b k f 


(62) 


that is governed by the differential equation 

i k . = J f ( B I< F - h K f ) = B P F - 11 P F = < „ (63) 

where ( () is defined to he the error in power, or the difference* between the actual and required values. 
The LQR technique can he used to control the kinetic energy error by minimizing the cost function 

Ve = / (ActT T r p ~) dt (0-4) 

Jo 

where A is a weighting parameter on the kinetic energy error of the* flywheel system. This leads to 
a feedback controller in which the required power is adjusted by the amount -\f\ r*.. Thus, the 
commanded power is defined to be 


B P C F = B P F -V Ac*. (<> r >) 

and is used in place* of gq (43) for the pseudo-inverse steering law. Similarly, a commanded 

power B P ( 1, (z — 1,2,3) is obtained for each of the three flywheel pairs and used together with 
Eqs. (51) and (52) for the divided power steering law. The merits of this kinetic energy error 
feedback are illustrated presently. 


Control System Block Diagram 

The block diagram in Fig. 3 shows the LQR (controller), flywheel steering law. and kinetic energy 
error feedback arranged to form a complete control system in MATLAB/Simulink ^ . The differential 
equations governing the behavior of the plant model arc* placed in a Simulink'^ S-Function to be 
numerically integrated and provide the value of the state {;r } at the current simulation time. The 
current value of the state is fed to the controller S-Function, where it is math* liondimensional and 
then multiplied by tin* LQR-derived state feedback gain matrix to generate tin* control torques r 
and t. The controller also determines the kinetic energy error of the flywheel system, r* . needed to 
calculate the commanded power, B P ( h \ that is used together with the flywheel control torque r in 
the flywheel steering law to obtain the flywheel motor torque measure numbers in {;/} [see Eq. (42)]. 
These arc* supplied to the plant model, and the control loop cycle is repeated. (The subject of CMG 
steering is well understood: therefore, our simulations do not include a model of individual CMG 
gimhal motions.) 
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Figure 3: Control System Block Diagram 


TEA Seeking 

The Torque Equilibrium Attitude (TEA) of a spacecraft is defined as the orientation for which tin 1 
angular acceleration of D in N vanishes. The present ISS attitude and CMG moment uni control 
algorithm can keep tin* orientation in the neighborhood of a time-varying TEA almost indefinitely, 
without requiring any expenditure of propellant from the Reaction Control System. In what follows, 
the design of a control law that seeks TEA with CMGs and flywheels is described briefly, and 
then simulation results are presented that illustrate the performance and features of our control 
and steering laws. We first illustrate the consequences of failing to counteract damping, and then 
demonstrate the advantages of kinetic energy error feedback. 

Seeking a TEA requires regulating the states contained in the 18 x 1 column matrix 

{./•} = ! 0 $, 03, u t lr x H\ J h]dt* ] T (Ofi) 

Weighting matrices [Q j and [Z?j are constructed with the 1 associated values presented previously. The 
matrices [A]. [Q]. and [Z?] are dimensioned 18 x 18, 18 x 18, and 18 x 0. respectively, in accordance 
with the number of regulated state variables for this control sc heme. 

Simulation results for this control law, discussed in the remainder of this section, are obtained 
with the following initial values of the state variables. Angles describing the orientation of B in L at 
t = /() are 0\ (to) = 5° (pitch). 0>(O)) = 5° (vaw), and 0^(to) = 5° (roll). Angular speeds associated 
with B (with l uj b = 0) are •«,(*„) = -9.80 x 10“* rad/s, u 2 (V) = -1.12 x HP* rad/s. and 

//;db j) — 9.82 x 10“ i> rad/s. Rotor spin speeds are u.[(t 0 ) = Uo(to) = u$(to) — — 20, 000 ri)in. and 
//r,(0i) — iir(t 0 ) — a->(0)) = 20,000 rpm. Initial values of CMG momentum measure numbers arc 1 
Mb)) = lb 2 (to ) = M>u) =0. 


Damped Flywheel Rotors 

The detrimental effects of damping are illustrated by performing a simulation with C c ( ~ 10 _r> X- 
m-s. One might 1)(‘ temptc'd to neglect such a seemingly small effect, especially over the short term, 
but it is shown here 1 to be 1 troublesome if not dealt with over long periods. The performance 1 of the 
TEA-seeking control law without feedbac k of kinetic* energy error is recorded in Figs. 4 7. 

Figure 4 shows the time history of the attitude angles and the inertial angular velocity. The solid 
curve is used for 0\ (pitch), the dashed curve for 0> (yaw), and the dash-dot curve is used for 
(roll). The 1 average 1 value's of these 1 orientation angle's in the 1 steady state are referred to as average 
torque equilibrium attitude angle's, anil are approximately the 1 same as those shown in Red. [4j. -7.0° 
in pitch, — 1 . 2 ' in yaw, and — 0 . 2 ° in roll. The 1 lenve'r plot of Fig. 4 shows the 1 inertial angular velocity 
response 1 , with u j. to. and u A . shown with solid, dashed, and dash-dot curve's respectively. Although 
they arc 1 not shown, plots of the 1 magnitude's of h and H are 1 virtually identical, remain below 14.000 
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Figure (i: Flywheel Angular Speeds with Damped Flywheel Rotors 


Flywheel Motor Torque Flywheel Motor Torque 
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Figure 7: Flywheel Motor Torque with Damped Flywheel Rotors 
















N-m-s over the first orbit, and in the steady state remain well below 4.745 N-m-s, the rapacity of a 
single ISS CMG. 

The upper plot of Fig. 5 contains the actual power delivered to the spacecraft by the flywheel 
system as a function of time. The middle plot displays the error between the actual power and the 
required power due to the damping in the flywheel system, which leads to the large secular kinetic 
energy loss of more than 50,000 kJ after 10 orbits, shown in the lower plot. 

The kinetic energy error of the flywheel system is reflected in the angular speeds of the flywheels 
relative to D shown in Fig. 0. Speeds of the flywheel pair whose spin axes are parallel to b[ are 
contained in the* upper plots, with represented by the solid curve, and a- indicated by the dashed: 
speeds of the pairs whose spin axes are parallel to b 2 and b 2 are depicted in the middle and lower 
plots, respectively. The flywheel rotor speeds exhibit a secular decay. No significant difference 
appears between the behavior resulting from the pseudo-inverse and divided power steering laws, 
presented in the left and right columns, respectively. 

Time histories of flywheel motor torques, in the presence 1 of damping, are displayed in Fig. 7 
where the results in the left and right columns an 1 associated with the pseudo-inverse 1 and divided 
power steering laws, respectively. The effects of damping appear to be' negligible at first: however, 
after some 1 time it is apparent that damping causes the motor torque magnitude's to increase, with 
either steering law. Inspection of Eq. (M) indicates that the' secular decay in rotor speeds requires 
an increase in motor torque's in order to produce the required power. The' increase in motor torque 
magnitude's implies further power losses as more electrical power must be diverted to the motor- 
generators in order to meet the attitude control and power management requirements simultaneously. 


Counteracting Damping with Kinetic Energy Error Feedback 

To compensate for damping, the kinetic energy error feedback design is employed with a weighting 
parameter A of 1 s _ ". The* parameters of the previous simulation are used again, leading to the 
results reported in Figs. 8 10. 
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Figure 8: Flywheel Angular Speeds with Damped Flywheel Rotors and Kinetic Energy Error Feed- 
back 
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Figure 9: Flywheel Motor Torques with Damped Flywheel Rotors and Kinetic Energy Error Feed 
hack 
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Figure 10: Power. Power Error, and Kinetic Energy Error with Damped Flywheel Rotors and Kinetic 
Energy Error Feedback 












Figure 1 8 shows that kinetic energy feedback eliminates the secular decay of rotor speeds seen to 
result from flywheel rotor damping in Fig. (i. The flywheel motor torques are shown in Fig. 9. It is 
immediately clear that the magnitudes do not increase with time as they do in Fig. 7. The power 
error e p shown in the middle plot of Fig. 10 is quite small and leads to the kinetic energy error r*. 
displayed in the lower plot, which is small and periodic, in contrast to the secular decay obtained 
without energy feedback (A = 0 s -2 ). 

The kinetic energy error feedback method compensates for damping very effectively, as illustrated 
by Figs. 8 10; these 1 results are virtually the same as those obtained with C\i — 0 and no energy 
feedback, although space limitation prevent us from inc luding them separately. 


CONCLUSION 

General, nonlinear equations governing motion of a rigid spacecraft containing flywheels and CMGs 
arc 1 presented in vector-dyadic form. A set of twelve scalar equations is obtained by applying the* 
generic relationships to the special case of a complex gyrostat with three pairs of flywheels mounted 
in orthogonal directions. Existing literature contains equations for describing motion of a space- 
craft with CMGs; they follow from the generic ones under two reasonable assumptions. The exact 
equations for the complex gyrostat, and the approximate relationships associated with CMGs arc* 
combined to form approximate equations for a spacecraft with flywheels and CMGs. and subse- 
quently linearized and nondimensionalized in preparation for design of linear control laws. 

A control law has been designed for an Earth-pointing spacecraft, and numerical simulation shows 
that it performs well in controlling torque equilibrium attitude', the energy stored in counter-rotating 
flywheels, and angular momentum of the flywheels and CMGs. Two steering laws arc* developed 
for ensuring that attitude control and energy storage* requirements are met simultaneously by the* 
flywheels. The design of a method for feeding back error in the rotational kinetic energy of the 
flywheel rotors in order to eliminate problems caused by rotor damping is shown to be* effective 1 . 
Possible* topics for future* research include control in the face* of imperfect knowledge of inertia 
properties, and unequal distribution of the attitude control workload between the flywheels and 
CMGs. 

The promising results demonstrated here for control of power, momentum, and attitude of Earth- 
pointing spacecraft suggest that it would be worthwhile to examine the very important, class of 
inertiallv oriented spacecraft. 
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